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We study singularities in the I-V characteristics for sequential tunneling from res- 
onant localized levels (e.g. a quantum dot) into a one dimensional electron system 
described by a Hubbard model. Boundary conformal field theory together with the 
exact solution of the Hubbard model subject to boundary fields allows to compute 
the exponents describing the singularity arising when the energy of the local level 
is tuned through the Fermi energy of the wire as a function of electron density 
and magnetic field. For boundary potentials with bound states a sequence of such 
singularities can be observed. 
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1 Introduction 



Electronic correlations together with strong quantum fluctuations are known to determine the 
low temperature properties of quasi one dimensional conductors. Theoretical investigations 
using integrable lattice realizations of these Tomonaga-Luttinger liquids (TLL) together with 
numerical studies and field theory approaches such as bosonization have provided much of the 
insight into the peculiar properties of such systems. Experimental evidence for TLL behaviour 
on the other hand is still rare in spite of the tremendous progress in synthetization of quasi one 
dimensional materials or fabrication of nano-structures in which the transport of electrons is 
confined to a single one dimensional channel One reason for this is that the theoretical 

work on TLLs has concentrated on low energy bulk properties of perfect infinitely long system 
which are difficult to access experimentically. Recently, studies of the response of a TLL on 
local perturbations have become feasible due to the construction of integrable models with 
open boundaries and a better understanding of quantum field theories in the presence of a 
boundary. Local inhomogeneities may have a profound effect on the transport properties of one 
dimensional interacting electron gases, even leading to phases with vanishing transmission of a 
barrier [^|-|^. Finite size effects in the resulting open chains have been studied to understand 
the possible experimental consequences of TLL properties in such systems (see e.g. p|-|14|). 
Among possible consequences of local perturbations are Fermi edge singularities which may 
be observed in X-ray absorption amplitudes and — as will be discussed in this paper — in 
tunneling experiments [|T^-[T8|. In both cases the nature of the singularities is strongly affected 
by the properties of the TLL. 

In this paper we study tunneling from a resonant localized level into a TLL in the limit 
of low barrier conductance as observed in the current-voltage characteristic at zero bias. As 
a specific representation of the latter we choose the one-dimensional Hubbard model. In this 
lattice model we find — similar as in previous work on optical absorption processes — a rich 
spectrum of edge singularities due to the existence of bound states reflecting the separation of 



spin and charge in the TLL |^. An experimental realization of the tunneling processes under 
investigation might be a quantum dot (providing the localized level) coupled to a quantum wire. 
The reservoir supplying charges to fill the state in the quantum dot is left unspecified. The 
energy Ei of the local level can be tuned by varying a gate voltage of the dot. The only infiuence 
of the quantum dot on the Luttinger liquid considered below is the electrostatic interaction with 
its net charge (we consider the dot occupied with a single electron to be electrically neutral). 
This description is very similar to the model invoked to describe the X-ray edge singularity in 



1 



metallic systems |^T| and has previously been used to study tunneling from a resonant local 



level into two- and three-dimensional systems These considerations lead to the following 
Hamiltonian 

L-l L ^ 

+ - 66V(iVi,| + iVi,|) . (1.1) 

where 6^ (6) are canonical fermionic creation (annihilation) operators for a spin-| electron the 
localized state and creates an electron of spin a on site j of the one-dimensional chain. 
The chemical potential /i and magnetic field h = gfisH allow to control the filling factor and 
magnetization of the quantum wire. Upon variation of the gate voltage tunneling between 
the local level and the wire becomes possible if the energy Ei of the local level exceeds the 
Fermi energy. We restrict ourselves to the case where the barrier conductance is low. Hence, 
the transport is dominated by incoherent sequential tunneling processes and we can neglect 



Coulomb blockade effects and higher order processes such as "cotunneling" |23|. Within the 



"orthodox theory" |^ the current due to sequential tunneling is computed by application of 



the golden rule leading to 

I{Ei) oc Y,\iM,ib\0)\^S{En - ^0 - Ei) . (1.2) 

n 

Here \0) = b^O) denotes the ground state of the open Hubbard chain in the A^'e-particle sector 

with the local level occupied and hence vanishing boundary potential p. The sum in (|1.2|) 

extends over all eigenstates \n) of the chain in the (A^^e + l)-particle sector in the presence of 
the boundary potential p. Eq. ( |1.2| ) can be rewritten as a Fourier integral: 



oo 



I{Ei) oc Re / dt e'^' *(0|6^(t)ci,|(t)4 .(0)6(0)|0) (1.3) 



where E^ = E^ + iO. Near the threshold Ei ^ Eth the intensity exhibits a characteristic 
singularity: 

m) oc . (1.4) 

For non interacting electrons the exponent a can be expressed in terms of the phase shift at the 
Fermi surface . As in the case of the X-ray edge singularity one expects several thresholds if 
the electrostatic potential p is strong enough to form bound states in the TLL ||2^,^]. In this 
paper we want to study this problem for tunneling into a TLL where an additional dependence of 
the exponent a on the interaction parameters (i.e. electron density, magnetization and strength 
of the Hubbard interaction 4m) in ( |1.1D of the Luttinger liquid is to be expected from the results 



obtained for the related X-ray problem (see Refs. [|T9|, |2^, Q ) . In the following section we 
summarize the relevant properties of the model ( |1 . 1| ) obtained from its Bethe Ansatz solution. 
From this solution combined with results from boundary conformal field theory (BCFT) p8|-pO| 
we extract the spectrum of thresholds and the corresponding exponents a. 



2 Bethe Ansatz Solution of the model 

The Bethe Ansatz equations (BAE) determining the spectrum of 7i with empty local state (i.e. 
with boundary chemical potential]?) in the iVe-particle sector with magnetization M = ^Ne—N^ 
read 



e^'=^(^^+^)S.(A:,) = n j = -N,,...,K 

sm fc,- — A3 — lu 



5 (A ) ft ~ + ^ TT ^a- >^f3 + 2iu ^ ^ _^ ^ ^2 

j=\^c.-smkj -iu ^l-jl^ A„- A/3-2m' ■ ■ ■ ^ i 

where one should identify k_j = —kj and A_q, = — Aq,. The boundary phase shifts appearing in 
the BAE read 

Bm = ( f^) B.(X) = ^ . (2.2) 

\ 1 — pe*'^ / sm k — lu A — liu 

The energy of the eigenstate of Eq. ( |1.1| ) corresponding to a solution of the BAE is given by 

E = ^ (^/i - - - 2cos%j + hN^ . (2.3) 

In Refs. pT|-P^j the ground state and the low-lying excitations of this model where studied for 
small boundary fields. In the existence of boundary states for \p\ > 1 has been established. 
In the Bethe Ansatz solution these bound states manifest themselves as additional complex 
solutions for the charge and spin rapidities. In Fig. |l| the spectrum of bound states for u = 1 
is shown. Using standard procedures, the BAE for the ground state and low-lying excitations 
in the thermodynamic limit can be rewritten as linear integral equations for the densities Pc{k) 
and Ps(A) of real quasi-momenta kj and spin rapidities Aq,, respectively: 

with the kernel K given by 




+ K* \ (2.4) 



cos k a2u(sin k — A') 

K=\ I . (2.5) 

a2?i(A - sin k') -a4„(A - A') 



Here we have introduced ay{x) = ^ > f*9 denotes the convolution dyf{x—y)g{y) 

with boundaries A = k^^^ in the charge and A = X^^^ in the spin sector. These boundaries are 
functions of the external chemical potential fi and magnetic field h. Alternatively, in a canonical 
approach the values of k^^^ and A^^-* are fixed by the conditions 



fc(0) 



dkpc 



2 [iVe - CJ + 1 



X(0) 



dXps 



2 [N^ -Cs] + l 
L 



(2.6) 



-fc(0) ^ J-X(0) 

where Cc (Cg) denotes the number of complex k (A)-solutions present in the ground state 
The boundary phase shifts ( |2.2| ) and the presence of complex solutions to the BAE determines 
the driving terms and p° in ( ^.41) . Their explicit form can be found in Refs. [pO| , pT1-|33[| . 
Denoting the solutions of ( p.4|) without the constant contribution l/vr to the driving term by 
Pc and Ps we introduce shift angles 




1 + 2C, 



e' = -\L 

P 2 



X(0) 



A(0) 



(2.7) 



Following Woynarovich one can calculate the finite size spectrum of the model, reproducing 



the result of PT| : 

E = LCoo + foo + 

+ 



nv. 
~L 

TTV, 



1 



1 



24 2det'fZl 



(2.8) 



Here Lcqo and f^o denote the bulk and boundary energy, A^'^^ are non negative integers counting 
the number of particle hole excitations at the Fermi points and the Vc^s are the Fermi velocities 
of the massless charge and magnetic modes. ANc^s specify the quasi-particle content of the 
state, in a TLL these are holons/anti-holons in the charge sector and spinous in the magnetic 



sector of the theory. The dressed charge matrix Z [34-36 




(2.9) 



is defined in terms of the integral equation 



1 
1 



ecs(A') e..(A') 



(2.10) 



A 



3 Tunnel exponents 



Results from boundary conformal field theory allow to extract the exponent a in Eq. ( |1.4| ) from 
the finite size spectrum ( p.8| ) p8|-|30| : the Green's function of an operator O with dimension x 
on the complex half plane is given by: 

^A\0{T^)0\r,)\A) = (3.1) 

Conformal mapping of the half plane onto a strip of finite width L allows to extract the 
scaling dimension of the boundary changing operator O from the finite size spectrum ( p^.8| ) by 
taking differences of the energy of the system's ground state \A) in the A^g-particle sector 
without boundary potential and the energy of the lowest excitation \B, n) in the [N^ + 1)- 
particle sector with boundary potential p and non vanishing form factor |(i?,'n,|(9"''|y4)p (see 
Refs. [^,^,^,^ for details on the CFT approach to the asymptotics of correlation functions 
in the Hubbard model). In the present problem E\ is obtained from ( p.8|) by choosing /S.N^ = 6q 
while for E'^ we have to choose ANj^ = 1 + and AN^ = 6q corresponding to an extra spin-f 
electron created by (9^. Writing the 0{L~^) terms in the energy difference as (tt/L) {vcXc + VgXs) 
the corresponding edge exponent in Eq. ( |1.4|) is given as 

a = 1 - 2{xc + Xs) . (3.2) 

We now want to study the exponents at the several possible thresholds. To gain some 
more insight into the role of the boundary states we begin with a discussion of noninteracting 
fermions. 

3.1 Ferromagnetic case 

For sufficiently large magnetic field the electrons are polarized ferromagnetically, hence an 
explicit expression is available for the wave function in terms of a Slater determinant of single- 
particle states. Considering \p\ < 1 these are plane waves corresponding to real wave numbers 
k exist and we expect a single threshold. The corresponding edge exponent is a function of p 
and the density of electrons rie = N^/L (see also Ref. [0): 



a'- = 1 - 2x; = 1 - (^1 + ^arctan ^) + " ^^'^^ 

For IpI > 1 the boundary potential can bind an electron which leads to the existence of two 
thresholds, depending on whether the bound state is occupied in the final state. For an empty 



bound state we obtain the result of Eq. ( |373| ) for the edge exponent. A different exponent 
is found if the bound state occupied: 

= 1 - = 1 - (^^ arctan tan + y ) ^ ■ (3-4) 

These predictions can be checked by studying the finite-size behaviour of the form factors. 
From the conformal mapping mentioned above one expects |2y,^| 

mc\M - (^y and mc[m ~ (^^y . (3.5) 

where \p) and \p) denote the ground state and the lowest state with empty bound state in the 
(A^e + l)-particle sector with boundary potential p > 1. Note that exponent vanishes in the 
limit p 0. This coincides with exact result, I{Ei) oc -w 1 — j-6 {Ei + 2 cos(7rne)), which does 
not exhibit a singularity. 

Eq. ( |1.2| ) is now evaluated numerically. To avoid the use of explicit representations of the 
(5-function we use the integral 

J{E) = I dE, I{E,) . (3.6) 

J ~oo 

Typical numerical results for J{E) are shown in Figs. |^ and ^. To make the numerical analysis 
of ( |1.2| ) feasible we had to restrict the sum to the most important states. The error of 
this approximation has been estimated using the sum rule J{E Emax) = 1 — {^1)1 which 
is satisfied to >99% in all cases. After smoothing of J{E) and numerical differentiation it is 
possible to obtain I{Ei). The result is shown in Fig. ^ for several potentials. For positive p the 
exponent is always positive at the absolute threshold and negative at the second one for p > 1. 
The situation changes completely for p < where the ground state is always parametrized by 
real wave numbers k giving a negative exponent a at the absolute threshold. Occupation of 
the anti-bound state corresponding to a complex k leads to a positive exponent. 

For the X-ray edge problem one can show that the functional dependence of J{E) is nearly 



unchanged by increasing the system size I^S],^. This allows to extract quantitative infor- 
mations of rather small systems {L = 80 in the present case) by fitting of J{E) to the trial 
function: 

f{E) = aiE + br . (3.7) 

The resulting exponent c can be compared to the BCFT-results. Using the 600 lowest states 
contributing to the sum (|1.2| ) we obtain a good agreement with the CFT-results (see Fig. |[) 
for p > 1 and electron densities He ~ 0.4. For larger densities a bigger discrepancy between the 
numerical results and the BCFT predictions is found due to stronger finite size effects. 



The second threshold due to the presence of a bound state is most pronounced for p < —1. 
Here the jump of J{E) characteristic of a positive edge exponent occurs at the second threshold 
(see Fig. H). While the BCFT result for the edge exponent at the absolute threshold is aabs = 

= —0.33 the fit to the numerical data gives aabs = —0.46 — this indicates that a genuine 
singularity is strongly affected by finite size effects. On the other hand the numerical value for 
the positive exponent at the threshold corresponding to occupied anti-bound state, a'^ = 0.988, 
is in very good agreement with the CFT-result a'^ = 0.977. 

3.2 Magnetic field dependence of edge exponents 

For vanishing magnetic field one has A*^°^ = oo allowing to solve the spin part of the integral 
equations by Fourier transformation. As a consequence the dressed charge matrix Z ( p.9|) is 
function of a single variable ^ = ^{k^^^) |3^, which is defined by the following integral equation: 

^{k) = 1+1 dk' cos k' G(sin k - sin k')i{k') , (3.8) 
J-m 

where G{\) = ^Re (l + i^) - * (| + i^)] (^(a;) is the digamma function). Further- 
more, one finds 9^ = \0p yielding 

aabs = 1 - 2(x, + X,) = ^ - ^ (1 + - ^^)' (3.9) 

for the exponent at the absolute threshold. In Fig. |^ we present regions where the exponent 
aabs is positive as a function of electronic density rif. and boundary potential p together with 
the density dependence of the exponent for some fixed values of p. In a finite magnetic field h 
the bulk state of the Hubbard model is ferromagnetic below a critical particle density ric- This 
density can be calculated from 



h = — I dk' cosfc' ^ ^ . ■ (3-10) 



2u f^"^" , . .cos A;' — cos(7rnc 
— / dfc' cosfc' . ; 

For non vanishing magnetic field we will only consider electron densities above Uc- For h > 
the exponent is given by: 

, [(1 + - e;) z,, - (gg - ej) z,,] ' + [(gg - g-) z,, - (i + - g-) j ' 

aabs = I . (3.11) 

In Fig. 1^ the magnetic field dependence of this exponent is shown for several values boundary 



potentials p. Note the characteristic change of this curve near p = pi = u + vtP+T where 
a low lying excited bound state for a charge and a spinon (corresponding to a complex quasi 
momentum k and a complex spin rapidity A in the set of roots of ( |2.1| ) [ |20| | ) appears, see Fig. 



For h approaching the saturation field ( p.lO| ) the exponents can be calculated explicitly 
giving 



m 



Oiabsine ricih)) > for < p < pi and aabsirie ndh)) < for p > pi . 

The difference of the limiting values aXp = pi is /S.aabs = 1- The crossover due to this behaviour 
is clearly seen in Fig. |^. As before additional edge singularities arise as a consequence of bound 
states in the boundary potential. The corresponding exponents can be computed from the 
Bethe Ansatz equations as above (for details see Ref. I^O])- Each of the boundary states seen 
in Fig. Ill gives rise to a singularity ( |1.4| ) in the I-V curve. 



4 Summary and Conclusion 

We have studied the I-V characteristics for tunneling from a resonant localized level into a one 
dimensional interacting electron gas described in terms of a Hubbard model (|1 . 1| ) . Compared 
to tunneling into a higher dimensional system one finds a rich structure of thresholds due 
to the presence of various bound states in the many particle spectrum, each of them leading 
to a possible resonance for the tunneling charge. Their appearance may be understood as a 
consequence of the separation of charge and spin degrees of freedom of the electrons allowing 
to generate a current either by the electrons decaying into their holon and spinon constituents 
in the wire or by holons alone leaving the spinon bound by the electrostatic potential of the 
tunnel contact. Furthermore, the electronic correlations within the quantum wire strongly 
affect the nature of the singularities, i.e. the exponents a in ( |1.4| ). While we have considered 



correlation functions strictly at zero temperature, the above singularities are still observable at 
small finite T: for sufficiently long wires the current at the threshold depends on temperature 
as I{Eth) oc T^". Hence, the dependence of a on parameters such as filling factor of the wire or 
the external magnetic field should be accessible experimentally thereby allowing to determine 
the properties of the potential due to a vacancy of the localized level. 
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Figure 1: Spectrum of boundary bound states of (|1 . 1| ) for m = 1 with chemical potential 
H = 0.5 and magnetic field h = 0.3, for (a) p < and (b) p > 0. The thresholds pi for binding 
a charge and spin and p2 for binding of a singlet pair of electrons are indicated. 
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Figure 2: Numerical results for J{E) for a system of L = 80 sites with a density rif, = 0.2 of 
spin-| electrons and p = 3. The dashed hne is the fit to Eq. (|3.7|) . The jump of J{E) at the 
threshold is characteristic to a singularity with a positive exponent a. Its height is given by 
the matrix element (p|cj JO) which vanishes in the thermodynamic limit according to Eq. ( p.5|) . 
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Figure 3: Same as Fig. but for p = —3. 
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Figure 4: Numerical results for I{Ei) for systems of size L = 80 with 16 electrons for different 
p < (a) and p > (b). The singularities are suppressed due to the numerical differentiation. 
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Figure 5: Comparison of the fit ( |3.7D to the numerical data with the BCFT results (solid 
lines) for three different electron densities. 
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Figure 6: Region with positive absolute exponent Oabs resulting in a singularity in the I-V 
characteristics for vanishing magnetic field for fixed u = 1. The insets show aabs for fixed values 
of p as a function of the electron density Ue- 
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Figure 7: Magnetic field dependence of the edge exponent Uabs ior u — 1 and fixed chemical 
potential /j, = —0.2 for several values of the boundary potential p above and below the threshold 
Pi = 1 + \/2 for creation of the holon/spinon bound state. 
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